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Abstract 

In a physical renormalization scheme, gauge couplings are defined directly 
in terms of physical observables. Such effective charges are analytic functions 
of physical scales, and thus mass thresholds are treated with their correct ana- 
lytic dependence. In particular, particles will contribute to physical predictions 
even at energies below their threshold. This is in contrast to unphysical renor- 
malization schemes such as MS where mass thresholds are treated as step 
functions. In this paper we analyze supersymmetric grand unification in the 
context of physical renormalization schemes and find a number of qualitative 
differences and improvements in precision over conventional approaches. The 
effective charge formalism presented here provides a template for calculating 
all mass threshold effects for any given grand unified theory. These new thresh- 
old corrections may be important in making the measured values of the gauge 
couplings consistent with unification. 



1 Introduction 



Precision measurements of the gauge couplings and their possible unification provide 
some of the few windows to the Planck scale. It is thus important to have a firm 
grasp of the theoretical ambiguities involved. This paper attempts to address some 
of these ambiguities. 

In a physical renormalization scheme, gauge couplings are defined directly in terms 
of physical observables. Such effective charges are analytic functions of physical scales, 
and thus the thresholds associated with heavy particles are treated with their correct 
analytic dependence. This is in contrast to unphysical renormalization schemes such 
as the MS scheme where mass thresholds are treated as step functions. In this 
paper we will analyze supersymmetric grand unification in the context of physical 
renormalization schemes with the goal of systematizing the effects of light and heavy 
mass thresholds and improving the precision of tests of unification compared with 
conventional approaches. 

In section 2, we motivate physical renormalization schemes with a simple example 
and then present the notation and results used throughout the paper. In section 3, 
we look more carefully at the problem of decoupling heavy particles and the errors 
induced by unphysical schemes. In section 4, we discuss the canonical self-energy-like 
effective charges for the minimal symmetric standard model (MSSM). These effective 
couplings run smoothly over spacelike momenta, have non-analytic behavior only at 
the expected physical thresholds for timelike momenta, and more directly measure 
the strengths of the forces than the charges of unphysical schemes. The extraction of 
effective charges from low energy data is considered. We identify an important modifi- 
cation of the electromagnetic coupling aq^D^Mz) due to the proper inclusion of virtual 
loops, thus resulting in a 4(T change in its numerical value. Similar modifications 
are found for the weak mixing angle. As seen in section 5, these effective charges pro- 
vide a more natural and physical framework for examining gauge coupling unification. 
In section 5.1, we demonstrate the invalidity of neglecting heavy threshold corrections 
in analyzing grand unified models. The more rigorous treatment of light thresholds 
in physical schemes gives rise to new corrections, but these are numerically small for 
most sparticle spectra. The treatment of heavy thresholds with various unification 
boundary conditions is discussed in section 5.2. In the simplest scenario, we find that 
the gauge couplings should unify at asymptotically large energies and the only heavy 
threshold corrections are logarithms of heavy mass ratios, corrections which can be 
obtained in unphysical schemes. An effective unification scale, defined in section 5.3 
as the scale where quantum gravity corrections produce non-negligible splittings be- 
tween the gauge couplings, is found to be roughly 10^'' — lO^^GeV, depending on the 
specific GUT model used. Section 5.4 considers more general unification boundary 
conditions with finite unification scale. The resulting heavy threshold corrections are 
given in Eq. ()44j) . This result combined with the results of section 5.2 may be used to 
determine the experimental consistency of any given GUT model. This is the main 
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result of the paper. An appendix discusses the details of constructing the effective 
charges. Throughout our analysis, we will find several attractive theoretical features 
of the supersymmetric regulator, dimensional reduction, which makes it the preferred 
regulator for physical effective charge schemes, even without supersymmetry. 

There have been several previous works on threshold effects in grand unification. 
In the first such study [T], which appeared just after the discovery of the grand 
unification, D.A. Ross uses form factors to define beta functions which are valid 
over all energy scales, including near mass thresholds. The coupling constants run 
smoothly over all momenta, and nontrivial threshold corrections are found for grand 
unification. Despite this early significant work, most subsequent work on GUTs have 
ignored these threshold effects, perhaps due to the complexity of the Ross approach. 

An exception from the late 1980's is the work of Kennedy and Lynn [3 , who 
defined electroweak effective charges similar to the pinch technique charges used in 
this paper. 

In several papers [2] by Kreuzer, Kummer, and Rebhan, the authors compared the 
Vilkovisky-DeWitt effective action (VDEA), the mass-shell momentum subtraction 
scheme(MMOM), and Weinberg's effective gauge theory (EGT). They wrote down ex- 
plicit formula for the running charges which include analytic threshold behavior for all 
particles. In calculating predictions from grand unification, they assume asymptotic 
unification at energies much larger than heavy particles, so that the only threshold 
corrections from heavy particles come from finite constants which are independent 
of energy scale or masses. We find similar results in section 5.2. Furthermore, we 
include the possibility of a finite unification scale in section 5.4, which leads to more 
complicated corrections. 

In [4j 5^ , the authors include the effects of light supersymmetric scalar and fermion 
thresholds, although heavy thresholds and gauge bosons virtual effects are not treated. 
In jS], the authors include both light and heavy threshold corrections, although the 
treatment of gauge bosons is not adequate. In Refs. [71 |H], the authors come to 
several conclusions similar to ours. However, their definition leads to gauge parameter 
dependent effective couplings. 

2 Physical Renormalization Schemes and Effective 
Charges 

In order to motivate the re-analysis of supersymmetric unification given in this paper, 
we will first discuss some general properties of renormalization schemes in the presence 
of massive fields and determine a criterion for consistent physical renormalization 
schemes. These criteria will not be satisfied by the schemes conventionally used in 
unification (and most perturbative calculations), MS and DR, which have persistent 
threshold and matching errors. Heuristically, these errors can be understood by noting 
that such schemes implicitly integrate out all masses heavier than the physical energy 
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scale until they are crossed, and then they are "clicked" on with a step function. 
Of course, integrating out heavy fields is only valid for energies well below their 
masses. This procedure is problematic since it does not correctly incorporate the finite 
probability that the uncertainty principle gives for a particle to be pair produced below 
threshold. Effective chargep^ schemes, derived from physical observables, naturally 
avoid such errors and are formally consistent. 

2.1 A Simple Example 

For the purpose of elucidating the benefits of physical renormalization schemes, we 
will give a simple toy example using QED with three fermions, e, /x, and r. Consider 
the amplitude for the process e~/i~^e~/i~. This can be written as a dressed skeleton 
expansion, i.e. the dressed tree level graph plus the dressed box diagram plus the 
dressed double box, etc.. The tree level diagram, dressed to all orders in perturbation 
theory, is equal to the tree level diagram with one modification : the QED coupling 
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a = 1^ is replaced by the Gell-Man-Low-Dyson effective charge 

"^^'^ " i + n^^(g2)_n^^(o)' 

Hence, from measurements of the cross section, one can measure the effective charge 
at two different scales, a{QD and a{Qf). Suppose the value of the electron charge 
is not known, and we are trying to test the predictions of QED. The way to proceed 
is to use one measurement, say at the low scale Qi, as an input to determine e. 
Now the prediction at the high scale Qh is well defined, and represents a test of 
the theory. More directly, we could just write a{QD in terms of a{Qf), leading 
to the same prediction. Since the cross section o^e-^i-^e-^i-{Q^) is proportional to 
{a{Q'^)Y, we are clearly relating one observable to another. The procedure just 
outlined is simply an on-shell renormalization scheme if Qi = 0. More generally, 
we will refer to such a scheme as an effective charge scheme, since we are writing 
a given observable, here just ae-^-^e-^l-{Q\) (or a{Q\)), in terms of an effective 
charge, a{Qf), defined from a measurement of the cross section at the scale Qi. One 
could equally well write any observable in terms of this effective charge. Note that 
this approach to renormalization works for arbitrary scales, even if the low scale lies 
below some threshold, say Qi < m^, while Qh > m^-. Decoupling and the smooth 
"turning on" of the r are manifest. 

Now we will compare with the results obtained by using the conventional imple- 
mentation of MS, which is as follows. First, the cross section is calculated at Qi 
using the rules of MS, which allows only the electrons and muons to propagate in 
loops, since Qi < rrir. Comparing the observed cross section to this result will fix 
the value of the MS coupling for two flavors, a2{,Qi)- To predict the result of the 
same experiment at scale Qh > rn-r, we need to evolve 0.2 to the tau threshold using 
the two flavor beta function, match with a three flavor coupling, ds, through the 
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relation 62(^1,-) = 0:3 (mr), and then evolve 0:3(771^) to Qh using the three flavor beta 
function. We will now have a prediction for o^e-n-^e-^l-{Qh) 0*^ ('^(QD)^- might 
expect, from the general principle of RG invariance of physical predictions that this 
result should be the same as the prediction derived using the physical effective charge 
scheme above. However, there is a discrepancy arising from the incorrect treatment 
of the threshold effects in MS. A detailed discussion of this problem will be given 
in section 3. In any case, the result can be obtained by straightforwardly applying 
the procedure outlined above. One flnds that the ratio of the cross section derived 
using MS with the cross section derived using effective charges, to flrst order in 
perturbation theory, is given by 



"(Ql) 3,1 

where Lr is a logarithm-like function (the high energy limit is a logarithm) given by 

LriQ/m) = I (ix6x(l — x) log ( 1 + — x) ) + 5/3 

JO \ j 

= (/3tanh-^(/3-i)-l)(3-/52) + 2, (3) 



where (3 = It satisfles the property Lr{0) = 5/3, so that there is no 

discrepancy when the low reference scale Qi is much lower than the tau mass threshold. 
This reflects the important, but often overlooked, fact that unphysical schemes, such 
as MS, are formally consistent only in desert regions where particle masses can be 
neglected. The error is plotted in Fig.(P). Notice that in this example there is an error 
only for Qi < nir. However, in the more general case of multiple flavor thresholds, 
there will be errors from both high and low scales. Similar discrepancies will be found 
in our analysis of grand uniflcation. 



2.2 General Properties of Effective Charges and Physical 
Renormalization Schemes 

Effective charges [HI may be deflned for any perturbatively calculable observable 

0{Q) = A'^ + a^a''{Q) (4) 

by absorbing all of the radiative corrections into the effective charge a'^ . To one- loop 
order using dimensional regularization (DREG) or dimensional reduction (DRED), it 
is straightforward to show that any unrenormalized effective charge may be parame- 
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Figure 1: The error in the MS based prediction for the scattering cross section, 
100% X (^1^ — l) , plotted against the reference subtraction scale Qi for the choice 
a{Qi) ^ 0.1. 



terized as^ 

2 2 

a'^iQ) = ao-^Y.^{L,{Q/m,)-r^'^{Q/m,)-Cuv + \og'^)+■■■, (5) 

where the sum is over all particles p in the fundamental theory which contribute to 
the running of the effective charge. In the QED example above, the sum proceeds 
over e,fj,,T and L^, = = L^-, and the function rj^{Q/mp) = 5/3 is a constant for 
the simple observable O = ^ae-f^-^e-fi-- In Eq-®, Cc/y = ^ — 1e + log47r is the 
divergence and associated constants, fi is the regularization scale, ao is the bare gauge 
coupling, and Pp is the contribution of each particle to the one-loop beta function 
coefficient. The Lp{Q/mp) are logarithmic-like functions which are characteristic of 
the spin of each particle, and are given exactly in Eq. (j23|) . They may be approximated 
to within a few percent^ by 

Lp{Q/m) ^ log (e^" + ^) (6) 

"'^This follows from considering the high energy limit and requiring renormalizability. Note also 
that our parameterization can be easily extended to effective charges which have particles with 
different masses running together in the loops, and the results are similar. In any case, we will not 
have use for such charges in this paper. 

^To be precise, the approximations reproduce the exact functions Lq, Li^2j Li (the subscripts 
refer to the spin of the massive field) with maximum error of 3.5%, 0.8%, and 2.2%, respectively, 
over the entire range of Q. 
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and have the hmits 



Lp{Q/m) ^tT' log^, Lp{Q/m) 



(7) 



where the constants rip have values given in the table below. We will see that these 
constants are of central importance in physical renormalization schemes. These log- 
like functions characterize the self-energy-like effect of each particle, including the 
finite spread of the wavefunctions near thresholds due to the uncertainty principle, 
and may be calculated in several different ways, as will be discussed in section 4 and 
the Appendix. Figure 2 shows the Lp functions for spacelike momenta. 



The Logarithm— like Functions 
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Figure 2: The logarithm-like functions for massive particles of spin 0, 1/2, and 1 are 
denoted by Lg^Lj, and Lw, respectively. 



Table I. 





scalars 


fermions 


massive gauge bosons 


Vp 


8/3 


5/3 


40/21 



The functions r]p{Q/mp) are characteristic of each observable, with a nontrivial func- 
tional form indicating deviations from self-energy like behavior. For a general ob- 
servable O, the function rj^iQ/irip) is nontrivial. We will show in section 4 that the 
constants rjp correspond to a particularly simple and canonical observable, called the 
pinch-technique (PT) self-energy. 

The effective coupling renormalized in the most general scheme R is 

= c.R{Qo)-^^^^§^Y.^(Lp{Q/mp)-Lp{Q,/mp) 

41:71 p 

- V?{Q/mp) + ii^{Qo/mp)), (8) 
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where the functions r]p{Qo/mp) contain all of the information about the scheme. Here 
R can be any mathematical scheme for defining the couplings. In the case of MS, 
we have rj^^ {Qo/nip) = Lp{Qo/mp) — log (Qo/^^p) so that only logarithms of the 
renormalization scale, Qo, are subtracted.^ 

It is straightforward to relate observables to each other: 

41:71 p 

- V?KQi/mp)+V?'iQ2/mp)). (9) 

This satisfies the transitivity property of the physical renormalization group. As be- 
fore, the sum over p runs over all particles in the fundamental theory which contribute 
to the effective charges. 

For consistency, very massive particles must decouple properly and must not con- 
tribute to physical predictions. Taking the m^— >cxd limit in Eq.(jH)) and Eq.(jni) yields 
a fundamental requirement of renormalization schemes and observables: 

<(o) = <H0) = <H0) ^^=^^ Vp- (10) 

This consistency requirement holds for all schemes R, observables Oi,02, and for each 
massive particle p. These are universal constants for each spin and, when DRED is 
used, are equal to the rjp given in Table I above, as can be verified through explicit 
calculations. If, for example, DREG was used instead of DRED, the last equality 
of Eq. p()|l would hold only for fermions and scalars, but for spin 1 fields there is an 
additional constant.^ Renormalization schemes that satisfy (the first two equalities 
of) Eq. lfTm) will henceforth be referred to as physical renormalization schemes, and 
those that do not will be called unphysical renormalization schemes, for reasons that 
will become clear. 

The above discussion implies a unique decoupling limit (Q/m^O) for observables. 
It is interesting that there is also a restriction on the high energy behavior (Q/m— >oo), 
which holds only for supersymmetric theories and takes the form of a sum rule. It is 
given by 

Eps5/?p(gX(oo) ^ ^^^^ 

Y.p(is(^p{G) 

where K'-' is a constant that depends only on the observable, not on the gauge group 
G or the super mult iplet S. The rjp{Q /nip) are calculated using DRED, otherwise the 
sum rule is true only for differences ?7^^(cx)) — rjp'^ {oo) between observables. Further, 

■^The term in parentheses in Eq.JSJ becomes Lp{Q/nip) — log ((3o/™p) ~ Vp{Q/''^p)- Note that 
in most calculations the first term is taken to be a logarithm and mass corrections are systemat- 
ically added, in order to approximate the full threshold dependence of Lp{Q/mp). However, the 
log {Qo/rrip) term does not have the correct threshold dependence, as we will be discussing. 

"'This is explained below Ea. lfT^ and below Ea. (|5S|l . 
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the result holds for any number of supersymmetries, which may be broken or unbroken 
at low energies. This can be proven inductively given the result for = 1. It is 
easy to check using Table I above and the corresponding result for massless gauge 
bosons given below Eq.jini) that K^^ = 2. The sum rule just expresses the fact that 
there is no resolution within a supermultiplet at high energies, and is motivated from 
conformal invariance and physical renormalization scheme invariance. Such a sum 
rule may provide a powerful link between the contributions of various spin fields to 
any observable, particularly if a multi-loop or non-perturbative generalization was 
found. 

3 Un-Physical Renormalization Schemes and the 
Problem of Decoupling 

Notice that the physical renormalization scheme requirement Eq. ljlOj) is not met by 
MS, DR, their massive extensions, or similar schemes. It is well known that MS by 
itself does not constitute a complete scheme, rather one must truncate the sum over 
p to include only particles with masses less than the scale of the problem. For each 
region between thresholds a different scheme is implemented and one must translate 
between schemes when crossing thresholds. Hence, MS is really a set of schemes 
related to each other. We will call such a set an artificial decoupling scheme(ADS), 
and now discuss the most general case at one loop. This will give us an idea of the dis- 
crepancies one may expect in ADS's when compared to the physical renormalization 
scheme approach. 

Let S = {mi,m2,...} be the spectrum of massive particles of the fundamental 
theory ordered from lightest to heaviest, let 5*0 be the set of massless particles, and 
let S'^ = So(B{mi,m2, ...,miy} be some subset up to a given mass scale. For any 
given renormalization scheme R, let the A^th phase of R, denoted R^ , be the scheme 
used to renormalize observables at energy scales Q such that 



To use R , one simply renormalizes the contributions from particles p E S in the 
usual way dictated by scheme R, and then entirely neglects the contributions from all 
p E S — S'^ . This is a formal statement of the usual implementation of ADS's using 
step functions. 

For mi, ...,m]\r < Q,Q' < m^^i, the gauge coupling of the R^ scheme flows by 



mi,..., mN < Q < rriN+i. 



(12) 



a^iv(g) 



^ Pp[Lp{Q/mp) - Lp{Q'/mp) 




(13) 
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and the most general matching condition between schemes i?^ ^ and takes the 
form 

tt^ ( '] 

ajfiv-i(/i7v) = a^iv(/iAr) + —pNANifJ^N/fTlN), (14) 

4:71 

where Ap{fip/mp) is arbitrary now, but will be specified below by minimizing errors, 
and may depend only on the matching scale /Xp for each threshold rup, for reasons 
discussed below. 

For Qh > rriN+n and rrii, ...,mN < Qi < rriN+i (the h and / stand for heavy and 
light scales, respectively) one may relate observables by fiowing through n thresholds 
using the above formulas to obtain 



a \Qh> = « [Qi> - 



E Pp{Lp{Qh/mp) - Lp{Qi/mp) - rj^'iQh/ 



47r 

+ V?%Qi/mp))+ E (3p[Ap{fip/mp) - Lpifip/rrip) + ri^{fip/mp) 

p(zSN+n_sN 

+ Lp{Qi/mp) - r]p^{Qi/m. 



TUp) 



(15) 



Now let us compare this to the relation, Eq.(jni), obtained in the previous sec- 
tion for the tracking of two observables. For the case of a high scale desert region 
(Qft-CmAT+n+i), the first sum reduces to Eq.Q. However, when Qh ^ niN+n+i, there 
are errors of one loop order which are proportional to Lp{Qh/mp) — rip^{Qh/mp), and 
occur for each neglected threshold^ p&S — S'^'^^ such that Qh ^ rrip. These errors are 
naturally remedied in physical renormalization schemes, leading to heavy threshold 
corrections which will be of importance when grand unification is discussed later. 
There are also analogous light threshold corrections. The second sum in Eq. ()15j) con- 
tains extra terms which arise from the artificial decoupling and matching conditions 
at each threshold. These terms are also generally of order of the one-loop corrections 
and must cancel if the ADS is to be consistent (in the sense of giving reliable physi- 
cal predictions in relations between two observables). In general, these terms do not 
cancel, since the Qz-dependent terms cannot be cancelled by the choice of Ap, which 
depends only on the ratio ^p/nip. Suitably choosing Ap (see EqUHl below) leaves a 
term Lp{Qi/mp) —Tlp^iQiImp) for all p^S'^^"' — . Hence, we see that the high scale 
and low scale threshold corrections have exactly the same form; indeed, they have 
the same origin, namely the necessity of using an ADS, which arises from improper 
decoupling. 

In the case of a low scale desert region, Qi<^mp\/p&S^^'^ — , one finds that 
LpiQi/rrip) — i]p^{Qi/'mp)^0 (by Eqs. ()7ll0p ). and the light threshold errors are elim- 
inated through the choice 

Apifip/mp) = Lp{fip/mp) - ri^{fip/mp). (16) 



'similar errors occur if Qi'^rriN+n+i 
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In this case, notice that once the Ap are chosen suitably, the matching scale fip 
exactly cancels and there is no need to fix its value. However, for MS, this choice is 
equivalent to using Ap = and fip = nip, which is the matching scale typically used 
at one loop. One may object that even in non-desert regions the known anomalous 
matching threshold errors could be systematically subtracted off for each physical 
process considered (equivalently allowing Ap = Ap{Qi)). This is untenable, as it is 
the same as using a different coupling for each process, thus losing the remnants of 
universality left by ADS's (i.e. universality in each desert region). 

We have identified two potential problems in artificial decoupling schemes, which 
arise solely from the failure of the decoupling requirement given in Eq. , regardless 
of whether or not the scheme R has analytic threshold dependence.^ The low scale 
errors come from the matching conditions and are exhibited in the last two terms of 
Eq. (fT3j) . These are significant only when Qi^rrif^^i. The high scale errors occur when 
one is calculating an observable at an energy Q/i^m7v+n+i that is slightly less than 
masses that should contribute, but are cut off in an ADS. These two errors will give 
rise to light and heavy threshold corrections in unification, as discussed in section IV. 

In practice, both types of errors can often be eliminated through a "threshold shift- 
ing" procedure. This involves modifying the definition of i?^ by replacing Eq. ()12j) 
with mi, ...,mj\r < aQ < myv+i and making similar subsequent replacements, and by 
choosing a > 1 to be large enough so that the desired thresholds that are slightly 
above Qi or Qh are 'moved' below aQi, so that no matching need occur for those 
thresholds, since they are already implicitly included in the couplings. The limit of 
this procedure as a^oo leads to a formally consistent scheme where no matching or 
artificial decoupling is used, but due to the failure of decoupling, it requires inclusion 
of contributions from every particle in the (unknown) fundamental theory. This is 
the exact situation that caused us to introduce an ADS in the first place, since we 
did not want unknown and arbitrarily massive fields contributing to every physical 
observable (written in terms of the ADS scheme charge). It is true that such unknown 
contributions cancel in relations between observables, but the utility of the interme- 
diate ADS scheme is lost since it's coupling is ill-defined^. In many calculations in 
unphysical schemes such as MS, the "threshold shifting" approach may be used to 
yield physical predictions which are arbitrarily accurate by choosing a sufficiently 
large 'a'. However, the usefulness of this procedure depends on the details of the 
mass spectrum. There is no universal algorithm that applies to any field theory. The 

^Proper analytic threshold dependence may be defined by rj^{Q/m) going to a constant for both 
small and large Q/m. Consider the analytic extension of MS into the region of mass thresholds, 
which we call massive MS, or MMS (similar to that in TIF). This is defined by Tjp^^'^^ (Qa/mp) — 
so that the full logarithmic-like functions Lp are subtracted and trivially the conditions for smooth 
threshold dependence are satisfied. Nonetheless, MMS has matching errors, which result from the 
failure of Ea. l|l()|l and the subsequent need to construct an ADS from MMS. 

^We might as well always write observables in terms of other observables; this is precisely the 
philosophy of effective charge inspired physical renormalization schemes. 
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complicated nature of such artificial fixes to the decoupling problem are reflections of 
the unphysical nature of the schemes and couplings. See ^T] for another approach to 
fixing MS. 

Thus, we have shown that the MS and DR schemes suffer errors unless one is 
restricted to observables at energies Q which lie far between mass thresholds. In 
addition, complicated matching conditions must be applied when crossing thresholds 
to maintain consistency for such desert scenarios. In principle, these schemes are only 
valid for theories where all particles have zero or infinite mass, or if one knows the 
full field content of the underlying physical theory. 

4 The Canonical Physical Effective Charges for 
the MSSM 

The difficulties associated with unphysical schemes are circumvented in physical 
renormalization schemes (PRS) based upon effective charges. So far, we have given 
consistency conditions which are satisfied PRS's, but explicit examples have not been 
given. This is the topic we now take up. 

For any observable (9, we define an effective charge scheme i?'-', by 

7^f{Q/m,;)=r^'^{Q/m,X (17) 

which, after using Eq.®, is equivalent to 

a'^(Q) = a^o(Q), (18) 

thus motivating the terminology "effective charge". Here i?^ is the physical subset 
of all possible mathematical schemes. The canonical example for using an effective 
charge as a scheme is furnished in QED by the Gell-Mann-Low-Dyson charge, which 
can be measured directly from scattering experiments. The extension of this concept 
to non-abelian gauge theories is non-trivial^], due to the self interactions of the 
gauge bosons which make the usual self-energy gauge dependent. However, system- 
atically implementing the Ward identities of the theory allows one to project out the 
unique self-energy of each physical particle, resulting in a self-energy that is gauge 
independent, may be resummed to define an effective charge, and may be related via 
the optical theorem to appropriate cuts of differential cross sections. The algorithm 
for performing the calculation at the diagrammatic level is called the pinch-technique 
(PT) [ISI lEI IIS] [in] • The procedure is illustrated in Fig. ((3]) for the case of massless 
gauge theory, where momentum factors from internal gauge boson lines or vertices 
combine with gamma matrices to cancel internal fermion propagators, yielding a gluon 
self-energy-like graph. This is then added to the usual self-energy to yield the full PT 
self-energy. 
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The PT procedure is unambiguous at one loop and is merely an application of 
the Ward identities of the theory, which becomes more transparent in a dispersive 
derivation from physical cross sections (j{qq—^gg) (see ^HI for such a construction for 
the electroweak sector). The generalization of the pinch technique to higher loops 
has recently been investigated. ^Zj^HmHlEn] [2I][22]- In the work of Binosi and 
Papavassiliou [20] [U] |22| , the authors prove the consistency of the pinch technique to 
all orders in perturbation theory, making clear how to define the QCD and electroweak 
effective charges at higher orders. 



self-energy-like 
projection 




self-energy-like 
projection 



self-energy-like 
projection 




Figure 3: Pinch-technique for QCD at 1 loop. The unique gluonic self-energy- 
like projection of the vertex and box graphs yield terms which must be added to the 
conventional self-energy to get the PT effective charge. 

The PT charge, labelled by a, written in terms of the bare coupling may be 
calculated for arbitrary gauge theory, broken or unbroken, to be* 

a{Q) = ao - (3p{Lp{Q / rrip) - % - Cuv + Im) -\ , (19) 

where Vp^~^^^^{:^) — Vp the constants given in Table I for massive fields and 
®If particles of different mass propagate together in the loops, this formula is trivially modified. 
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rjg = 64/33 for massless spin 1 fields^. The fact that these rj^{Q/m) functions are 
constants is what makes the PT observable the most simple and natural choice for 
defining an effective charge scheme. More general physical effective charge schemes 
(see Eqs. fl8|9fl()j) ) have more complicated running due to the rj'^{Q/mp) terms. The 
calculation of Oi{Q) has been performed using dimensional reduction (DRED), rather 
than dimensional regularization (DREG). We will let PT stand for the renormaliza- 
tion scheme associated to the PT observable regularized using DRED. The results 
in DREG for s = 0, 1/2 fields are the same, but for spin one fields are given by the 
replacement rjw rj^y + 2/21, where the additional constant is due to the so-called e 
ghosts. Although the constant terms will cancel in relations between physical observ- 
ables, DRED is the more natural choice, even for non-supersymmetric theories. This 
is because decoupling is manifest (including for massive spin 1 fields) individually 
for each observable without relying on cancellation terms arising from relations to 
another observable. This is just a statement of the fact that f^P'^-^^^D _ /^^(Q)Yp_ 

In contrast, with DREG one has ri^~^^^^ ^ Lvi/(0). Also, the supersymmetric sum 
rule in Eq. ()ll|) becomes manifest in DRED. 

Using the above results, it is straightforward to write down the effective charges 
for the standard model through 



ai{Q 



2n 



5 a(Q2 



2\ 



3i-s2(g2) 

a{Q') 



asiQ') = as{Q'), (22) 

where the effective couplings a and are defined from PT self-energies Il^-y and II^z, 
respectively as is detailed in Appendix A. It is convenient to write 5i and 02 in 
terms of a and since the latter contain the contributions from the mass eigenstate 
fields. One could use Eq.lfT^ directly, although the Higgs sector requires care. 

Several subtleties should be addressed before the numerical values of the PT 
couplings are given. 



^Wc will use 'W' or '1' subscripts to denote massive spin 1 fields and a 'g' subscript for massless 
spin fields. The constants 64/33 and 40/21 are related straightforwardly. In general, for a massive 
gauge boson W in the representation R of group G that is being considered and representations R' 
in additional group factors G", we have 

Pw = ^C(R)d{R') - ^C{R)d{R') = ^-C{R)d{R') (20) 



and 

Vw 




c(mR'}Q-lc(m^(n'){l)]-i. (21) 
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An important difference between the physical effective charges and the unphysical 



MS couplings is a distinction between timelike and spacelike momenta. In conven- 
tional approaches, thresholds are treated in a step function approximation, and hence 
the running is always logarithmic. The analytic continuation from spacelike to time- 
like momenta is trivial, yielding in imaginary terms on the timelike side. Thus, the 
real parts of such couplings are the same modulo three loop (ivr)^ corrections. In 
contrast, the PT couplings on timelike and spacelike sides have considerable differ- 
ences at one-loop. To see this we need the exact expressions for the logarithmic-like 
functions of a particle of spin s, which can be written as 



Ls{Q/m) = 2 



/ACS _ ^2^ 

(/?tanh-i(/5-^)-l)U^^J+l 



(23) 



where S'^ = s{s + 1) is the total spin squared eigenvalue, /? = ylT^p-, and the 
momenta is spacelike (Q^ > 0). This formula is merely a compact way to write 
the results for massive spin 0, 1/2, and 1 fields, and has not been explicitly verified 
for higher spins. For example, Lg^i is calculated from the sum of the usual gauge 
boson self interaction loop, the ghost loops, the appropriate loops of Goldstone bosons 
that are eaten, and the pinched parts of the vertex and box graphs (see appendix A 
for details). In contrast, Ls=i/2 is simply the related to the usual fermion vacuum 
polarization graph, and Ls=o comes from the usual scalar contribution to the gauge 
boson self-energy. It is interesting that such a simple compact form is obtained, 
considering the seemingly different derivations of the three Ls functions. This may 
suggest a more efficient formulation of the perturbative dynamics of quantum fields 
that treats the various spins in a unified manner [SB]- Notice that 



lim LslQ/m) 



1 -3s(s + l) 
1-As(s + 1] 



(24) 



corresponding to the results of Table I. The analytic continuation of Eq. ()23p to time- 
like momenta below threshold, < = —Q^ < 4m^, is obtained by replacing 



/ 4?77'^ — 1 

where /3 = - 1, and tanh-^(/5-^)^ - itan~^(/? ). (25) 



Above threshold, > 4m^, one should replace 



tanh-\(3-^)^tanh-\(3) + where /5 = Jl - (26) 

Z V ? 

From these results it is clear that significant differences will arise between the spacelike 
and timelike couplings evaluated at scale M|, mainly due to the boson threshold 
asymmetry. 



Am? 
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As has been discussed, another distinction of effective couphngs is that they are 
automatically sensitive to light SUSY thresholds near Mz, since the Lg functions are 
not zero below threshold (on the spacelike side nor on the timelike side). The effects 
of light SUSY thresholds on the values of the couplings at the Z-pole will depend 
on the method of extraction from the data. The key question is whether or not the 
light sparticles are implicitly included in the measured values of the couplings at 
Mz- For a{Mz), which is extracted by running the precisely known fine structure 
constant from Q = to Mz, we should include corrections from virtual effects of 
sparticles (with model-dependent mass), in the self-energy term 11^^ (M^). However, 
these threshold corrections will cancel in any unification prediction, since then one is 
essentially running from Q = to Q = Mqut and the light SUSYs are either fully 
decoupled or fully turned on. For the strong and weak couplings we use data from the 
Z-pole, and thus no unknown sparticle thresholds must be accounted for since they 
are already implicitly contained in the measured values. When these couplings are 
run to the unification scale the induced light threshold corrections will not cancel. Of 
course, linear combinations of the electromagnetic and weak couplings (Eq.(j22I)) are 
used for unification, which complicates the matter further, since different methods 
of extraction are used for each. It would be unpleasant to quote a different value of 
a~^{Mz) for each different SUSY spectra considered. However, this approach has the 
advantage that the values of the couplings used are the values that one would directly 
measure in an experiment at Mz if a given sparticle spectrum were the correct one. 
For convenience, we will quote the QED coupling extracted assuming a fully decoupled 
SUSY. When calculating detailed unification predictions in given models, however, 
the appropriate terms will be included in the determination of a~^{Mz)- It should 
be emphasized that the above complications are only numerically significant for light 
sparticle spectra. 

The initial values may be extracted from experimental data and are given in the 
following table, where spacelike and timelike effective couplings are denoted with a 
and '— ', respectively. The MS and DR couplings are on the timelike side. 



Table II. 





MS 


DR 


PT+ 


PT_ 


a~\Mz) 


127.934(27) 


127.881(27) 


129.076(27) 


128.830(27) 


s\Mz) 


0.23114(20) 


0.23030(20) 


0.23130(20) 


0.22973(20) 


a^{Mz) 


0.118(4) 


0.119(4) 


0.140(5) 


0.140(5) 



See the appendix for detailed formulas for the effective couplings. 

Notice that the value of the PT+ electromagnetic inverse coupling, a~^{Mz) = 
129.076(27), does not correspond to the usual value of about 128.968(27). This dis- 
crepancy arises because a~^{Mz) includes the virtual effects of W~^W~ loops, whereas 
the usual construction of aqEniMz) entirely ignores the virtual effects of the massive 
gauge bosons. The proximate cause of this consistent oversight in the literature is the 
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difficulty in extracting a gauge invariant self-energy-like contribution to the running 
couplings for non-abelian theories, a problem which is resolved through the pinch 
technique, in particular, and more generally, in any effective charge scheme. Clearly, 
the present approach yields a coupling which more accurately reflects the strength of 
the electromagnetic force. Similar comments apply to the weak mixing angle. 

It should be emphasized that although we have chosen to discuss a particular 
physical renormalization scheme, it will be shown in the next section that all predic- 
tions associated with unification are PRS invariant, as they should be. However, a 
definite scheme must be chosen for explicit calculations, and the PT scheme is the 
simplest choice. As expected, we will find that PRS invariance does not extend to un- 
physical schemes such as MS or DR, because of errors associated with the incorrect 
treatment of light and heavy thresholds. 



5 Unification in Physical Renormalization Schemes 

Now we are ready to discuss unification. In section 5.1, we will consider only the light 
spectrum given by the standard model fields and their = 1 superpartners. This 
gives a model-independent starting point for discussing unification, and makes clear 
exactly what model dependent heavy threshold corrections are needed for consistency 
with the unification hypothesis. New light threshold corrections, in addition to the 
usual light mass corrections, are evident, although they are numerically important 
for only a small range of parameter space corresponding to light sparticles. In section 

5.2, asymptotic unification is introduced, leading to substantial qualitative changes in 
the usual picture of gauge unification. This particular choice of unification boundary 
conditions will lead to corrections from logarithms of superheavy mass ratios, just 
as would be obtained by implementing DR with the step function approximation. 
This sheds light on the nature of the approximation of the DR approach. In section 

5.3, an effective unification scale is derived that is considerable higher than the usual 
unification scale. In section 5.4, more general non-asymptotic boundary conditions are 
considered, and the new non-trivial thresholds corrections are found to be important. 

In performing the analysis, the exact analytic one-loop formulas discussed in sec- 
tion 4 will be used, as well as the leading two-loop corrections. The analytic mass 
dependent two-loop corrections are not known, but these can be estimated to be 
numerically small and well within the error bars, and hence can be neglected 

We will treat the SUSY spectrum as entirely arbitrary, rather than assume a 
particular model or theoretical bias. The advantage of this approach is that impor- 
tance of various spectra parameters becomes transparent, and irrelevant details can 
be ignored. 
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5.1 The (in) validity of Neglecting Heavy Thresholds 



In this subsection only, heavy thresholds will be entirely neglected. 

The usual test of unification is to predict a^i^Mz) contingent upon unification. 
Compared with the conventional DR framework, we expect to see improvements 
due to the correct treatment of light thresholds. To be precise, the corrections we 
are discussing are to the difference between the 03 (M^) prediction obtained from 
the following two methods: (a)using the PT+ scheme throughout, (b) using DR 
(with the artificial decoupling and theta function treatment of light thresholds) to 
predict a^i^Mz), which is then translated to a prediction for a^i^Mz). Both approaches 
capture the leading light threshold effects, which appear as logarithms of light masses. 
The additional corrections in the PT scheme are from what we will call analytic 
light threshold corrections, since they arise from correctly and smoothly interpolating 
between thresholds. These are largest when there are light supersymmetric partners 
near or below Mz- For most values of the sparticle masses, they fall inside the error 
bars. However, such corrections may become more important as the experimental 
values of the couplings are determined more precisely. The exact form of the new 
corrections will be shown explicitly in section 5.2, Eqs. (j31l35j) . 

Now let us compare the PT unification predictions with experiment. The predic- 
tions for the PT strong coupling, 03, (obtained through method (a)) are displayed 
in Fig.|3] against the SUSY scale and in Fig.|2] against the mass ratio of the gluino 
and wino. These are the two SUSY spectrum parameters to which the 03 prediction 
is most sensitive. 

Only light gluino scenarios with m~ ^ m~, are able to correctly predict the strong 
coupling for natural SUSY scales (less than about a TeV). However, it is generally ex- 
pected that the gluino is several times heavier than the wino for most realistic models 
of supersymmetry breaking and spectra. Hence, we reproduce the known result 
that, at two loops and neglecting heavy thresholds, gauge coupling unification fails 
by several standard deviations. Except for the light gluino escape route, this points 
to the need for large heavy threshold corrections if unification is to be achieved. 

5.2 Heavy Thresholds and Asymptotic Unification 

Henceforth, the complete heavy threshold behavior will be included in the running 
of the effective couplings. The form of the subsequent corrections will depend on the 
particular unification boundary conditions that are chosen, and the numerical values 
of the corrections will depend on the details of the GUT model. In this section we 
will choose the simplest boundary conditions, since it will reproduce known results. 
Later, more general cases will be considered. 

Generally, there are four parameters which specify the unification boundary con- 
ditions. These are the unification scale, Mu, and the values of the couplings at that 
scale, ai{Mu) for i = 1,2,3. For our purposes, we will always assume standard nor- 
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Figure 4: The error in the prediction for 03 (M^) is plotted against the typical SUSY 
mass scale, with different lines corresponding to values of the ratio of the gluino mass 
to the wino mass. The relative mass spectrum is roughly the same as most sparticle 
spectrum models, including supergravity models, with Mg setting the overall scale. 
The experimental standard deviation, s.d., is 0.0055 for the PT strong coupling. 



malizations and take the couplings to be equal at some scale. In this case, the only 
free parameter is M[/. The two distinct cases are for finite Mu and infinite Mjj. The 
so-called asymptotic unification considered in this section corresponds to the latter 
choice, namely Mu^oo and aj~^(Mt/) = a^^(M(/) = a^^{Mu). The asymptotic unifi- 
cation conditions would be appropriate if the standard model group Gsm is embedded 
in a simple Lie group G which is fully restored before gravitational or other string 
interactions become relevant, and neglecting any other exotic phenomena. Hence, 
this choice is somewhat simple and naive, but it is very instructive. 

We will find that asymptotic unification reproduces the same heavy threshold 
corrections which can be obtained by unphysical renormalization schemes (DR) with 
finite unification scale. The reason is that in taking Mjj^oo, one is essentially looking 
at an observable (the unification requirement) in a desert region, which, as we have 
seen, unphysical schemes are capable of treating without error. At first sight, it 
may seem strange that the infinite unification scale predictions of physical schemes 
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strong Coupling Predictions versus Gaugino masses 
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Figure 5: The error in the prediction for 53(71^^) is plotted against the ratio of the 
gluino mass to the wino mass, which is the sparticle spectrum parameter to which 
as{Mz) predictions are most sensitive. The spectrum is fully specified by the ratio 
and ^m~m~ = SOOGeV = Mgusy, where Msusy is the mass of all other sparticles. 

correspond to finite unification scale predictions of unphysical schemes. However, this 
is dictated by the nature of unphysical schemes where masses are turned on and off 
with a step function. 

The paradigmatic improvement over conventional methods is summarized in Fig. |B] , 
where asymptotic unification of the couplings occurs at very large energy. For demon- 
strative purposes, the parameters are chosen so that unification occurs. 

Now let us derive the analytic formulae for the unification predictions. We will dis- 
cuss the most general case of an = 1 supersymmetric Gsm = U{1)y®SU (2) l^SU {3)c 
embedded in a larger gauge group, G, using any physical RN scheme (all others are 
inconsistent), which we label by its associated observable, O. 

In general the running of the couplings can be expressed in the form 

c^oliQ) = <(Qo) + nf (Q, Qo) - 0^{Q, Qo), (27) 
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Asymptotic Unification 
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Figure 6: Asymptotic Unification. The solid lines are the analytic PT effective 
couplings, while the dashed lines are the DR couplings. For illustrative purposes, 
a-i,{Mz) has been chosen so that unification occurs at a finite scale for DR and 
asymptotically for the PT couplings. Here Msusy = 200GeV is the mass of all light 



M. 



SUSY 



2m- 



superpartners except the wino and gluino which have values ^m-- 
For illustrative purposes, we use SU{5). 

where the two-loop corrections^^ are contained in Oi{Q,Qo), and we have defined 
nf (Q,Qo) = 2^ E Pi'\L,iQ/mp) - L,{Qo/m,) - v?'{Q/m,) + r/f (Qo/mp)) , (28) 



p6G 



which contains all of the one loop corrections. Now we separate the sums over the 
light and heavy spectra, L = Gsm {Gsm means the standard model fields plus SUSY 
partners) and H = G — Gsm, fake Qo = Mz, and let Q = Mu be some energy much 
larger than the mass of all fields, including the heavy fields; i.e. Mu^irip \lp&L + H. 
The functions Iip can then be written as 

nf (M^-.oo, Mz) = PgIu - Af - 5f - (3flx - 5f,(oo) + 5f,(M^) - 5^_,(oo) (29) 
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see Appendix A for the details 



21 



and 



where /3g = Epec lu = ^ log §§Jx = ^ log 

sUQ) = T.^(^i'^v?{Q/^i), (30) 

^f^Ej^#'losg. (32) 

The exact 1-loop analytic hght threshold corrections are contained in Af, while the 
heavy threshold splittings are contained in , with some arbitrarily chosen heavy 
mass Mx which is conveniently taken to be the mass of heavy gauge bosons. 

It is useful to verify that predictions for Ix and a^lMz) are invariant under the 
choice of physical renormalization scheme. In performing the calculation, one must 
use the fact that the rjp functions do not depend on the gauge group or representation 
of p, only the spin. These are necessary (but not sufficient) conditions for the sum 
rule in Eq. ljllj) . This scheme equivalence does not extend to unphysical schemes such 
as DR, though the errors are quantifiable. 

Due to the physical renormalization scheme invariance, we may choose the simplest 
scheme, which is the PT scheme discussed earlier. Because the rj^^ functions are 
constants equal to rjp = Lp{0), the expressions for the unification predictions are 
simple and compact when written in terms of the PT charges 

From ai{Mu) = a2{Mu), the heavy gauge boson mass, Mx, is given by 



log {j^) + 1 a^\Mz) - a^\Mz) + A 



12 



27r (5: 



(33) 
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where A12 = Ai - A2, /5i2 = A - /32, and Ai = Af + 5f + 6*1. Notice that Mx can 
be determined explicitly only for the (unlikely) case of a degenerate heavy spectrum 
when 5f = 0, otherwise the expression is transcendental in Mx- In the degenerate 

case, the gauge boson mass Mx = J^/e is equal to the unification scale determined 
by entirely neglecting heavy thresholds(denoted by^), divided by e = 2.71828.... This 
result relies on use of the sum rule in Eq. ljlH) which gives rise to the l/27r term on 
the LHS of Eq. ()33|) . The generalization to arbitrary physical renormalization scheme 

is Mx = I\i(je~^° /"^ , where K'-^ is defined in Eq. (lll|) . Neglecting the light and 
heavy analytic threshold corrections, the gauge boson mass prediction is the same 
as the unification scale prediction in the DR scheme. Also, the 'unification' scale 

Aijj depends on the particular scheme, which makes sense since different schemes 
correspond to different observables. In contrast, Mx is scheme independent. 
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The strong coupling prediction is 

a^\Mz) = a^\Mz) + A3, + (^(5^\Mz)-a^\Mz) + Ai2), (34) 

Pl2 ^ ^ 

which differs from the prediction obtained by neglecting heavy thresholds by only the 
terms which reflect the heavy splitting. 

In order to explicitly compare with the DR approach, the artificial decoupling 
treatment of thresholds should be employed, as described in section 2. This involves 
using a step function through each light field / G L with mass greater than Mz, 
and through every superheavy field h & H. Then one must impose the unification 
condition that the three gauge couplings are equal at the maximum mass of heavy 
fields, Mx = max{'mh,h G H}. At energies above this maximum mass, the three 
couplings run identically according to the beta function for the unified group G; hence 
there is no arbitrariness in the choice of the unification scale. Next, the prediction 
for the DR strong coupling should be translated to the PT strong coupling. Doing 
this, one finds the exact same form of Eq. ljH^ . except that Af is replaced by 

mi>Mz V '"'I / 

Notice that there are only light threshold corrections beyond the theta function ap- 
proximation for particles of mass above Mz, since those below Mz are already im- 
plicitly accounted for. This formula is in agreement with Eqs. jlSpifij , since there is a 
residual error proportional to Lp{Mz/mp) — rjp for each crossed threshold. The anal- 
ogous corrections for the heavy thresholds do not arise in the asymptotic unification 
scenario, since we are essentially comparing observables at energy scales Mz ~ mi, 
which is of the same order of magnitude as the light thresholds, and Mu^rrih, which 
is much greater than all thresholds when asymptotic unification conditions are as- 
sumed. The latter scale is a "desert" scale, and so the step function method has 
no errors, giving the same result obtained above in the Sf^ .^^ For the more general 
unification conditions considered in the next subsection there will be additional heavy 
threshold corrections. 

Eq.(j21I) is a useful result, as it allows one to constrain the heavy spectrum, given 
a light SUSY spectrum. Up to two loop finite threshold corrections, which we have 
estimated to be small, and assuming that Eq. ()34|) will yield the experimental value 
of the strong coupling given some appropriate full GUT theory (i.e. assuming the 
asymptotic unification hypothesis is true), we can write 



e = a 



\Mzrr - as'iMzr^' ^ -Sg- ^Sii, (36) 



12 



""^""^It should be emphasized that this is only the case when the DR is correctly implemented by 
choosing the unification scale to be equal to the heaviest threshold in the theory. Different choices are 
sometimes made in the literature. For example, in ^2^, the authors advocate defining the unification 
scale to be the geometric mean of the heavy masses. 
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where a^, ^{Mz)^^'^ is the predicted value of the strong couphng obtained by neglecting 

heavy thresholds, as illustrated in Figs. j4p5j . We should emphasize the assumptions 
leading to this result. First, the standard normalizations of the couplings are assumed, 
so that Eq. ()37p does not hold for higher affine levels or non-standard hypercharge 
normalizations, as often occur in string models. Second, we assume that the gluino 
is somewhat heavier than the chargino, so that there are serious discrepancies, as in 
Fig. which must be explained by heavy threshold corrections. Finally, we are using 
the paradigm of asymptotic unification, wherein the full gauge group G in which the 
SM is embedded is restored before other Planck scale physics becomes relevant. With 
these assumptions, and noting that heavy thresholds were neglected in Figs. |4l5] . we 
find a typical value of 



«"-l"'-'-ii:B.log^. (37) 

where we have defined^^ 



= + (38) 

Values of can be compiled for the heavy representations any unified gauge group, 
and hence may be used with heavy mass ratios to exclude or provide evidence for a 
given GUT theory. 

To calculate B^, we first write B^ = P^^Bh, where = —1/3, —2/3, 11/3 for spin 
0, 1/2, 1 fields and the remaining group theory factor is Bh = |Ti(i?) — yT2(-R)+T3(_R) 
for a representation R. It is necessary to decompose all representations in terms of 
their U{1)y®SU{2)l0SU{3)c content. Here Ti(i?) = ^EpenV and Ti{R)5^^ = 
X]pgi?trj(t^t^), i = 2,3. For most grand unified theories of interest, all multiplets can 
be decomposed in terms of only eight different standard model multiplets (plus their 
conjugate representations, which have the same Bh, and a singlet which has Bh = 0), 
which are given in Table III along with the value of Bh- 



Table III. 





Bh{Ri) 




Bh{Ri) 


i?i = (3,2,1/6) 


-3/2 


i?5 = (1,1,1) 


3/7 


i?2 = (3, 1,-1/3) 


9/14 


i?6 = (8,l,0) 


3 


i?3 = (3,l,2/3) 


15/14 


i?7 = (1,3,0) 


-24/7 


i?4 = (l,2,l/2) 


-9/14 


i?8 = (3,2,5/6) 


3/14 



These same constants will also govern the corrections from analytic heavy thresh- 
old corrections that will be discussed later in section 5.4. Notice that, by definition, 

i^Note that = ^. 

P12 ' 
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the constants satisfy the constraint that the sum over all heavy multiplets vanishes, 

Y.Bh = Q, (39) 

which equivalently reflects the arbitrariness in the choice of which heavy mass scale 
Mx one chooses to be canonical (see Eq.(jnZj)), = 0. A similar relation also holds 
for any complete representation of the grand unified group. For example, the 24 of 
SU{b) decomposes into a singlet plus -Rg + -Ry + -Rs + -Rs- From the table, we have 
B,,{R^) + Bn{Ri) + 2Bf,{R^) = Q. 

As a simple example, let us explore the (unlikely) possibility wherein the only 
heavy field with significantly different mass than the heavy gauge boson mass Mx 
is the 5 dimensional Higgs supermultiplet in which the light Higgs doublets are 
embedded. The triplet components of the two Higgs supermultiplets contributes 
-2/5,0, -1 to pi,P2,p3, and hence Bh{3 + 3) = -9/7. Using Eq.(jSI|), this leads to 
M3 ^ Mx exp — , which is of order Mx/100. Such a large splitting is unnatural 
and difficult to accommodate in a theory. In general, "natural" splittings do not lead 
to values of the correct magnitude in SU{5). This is not terribly surprising, since 
minimal SUSY SU{5) is already known to be strongly disfavored. 

In general, the large discrepancies in Figs. j4l5] imply a large splitting in the 
heavy spectrum, which, in turn may imply a multistep unification scenario, e.g. 
SO {10)-^G224,-^G SM ■ The reason is that for the heavy fields to contribute to a^i^Mz), 
they must not only have a mass splitting compared to some reference heavy gauge 
boson, X, but also must have different first beta function coefficients since only the 
differences (3i^^ — (32^^ and (3i^^ — 13^'^^ appear in the corrections. 

Before moving to more general unification boundary conditions, we shall give 
a simple way to define an effective unification scale in the asymptotic unification 
scenario. 



5.3 Effective Unification Scale 

Because the couplings formally unify at infinite energy in the paradigm of asymptotic 
unification, there is no apparent unification scale. However, we suspect that in reality 
quantum gravitational fluctuations will affect the couplings as they approach the 
Planck energy. Hence, one can define an effective unification scale to be where the 
splittings between the gauge couplings are of the same order as those induced by 
gravitational effects. To be precise, define a dimensionless gravitational coupling 
which classically runs with energy as 

Gn{Q) = (40) 

where Mpi 1.22xlO^^Ge^. The leading gravitational corrections to the running 
gauge couplings ai{Q) will be proportional to GN{Q)o:i{Q). Hence, the effective 
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asymptotic unification scale, Meff, can be defined as the scale where the split- 
tings in the gauge couplings are of order the gravitational corrections, |aj(Me//) — 
aj{Meff)\ ^ b'^GN{Meff)au, or equivalently 

\a-\M,ff) - a-\M,fj)\ ^ b'GNiM,ff)au' ^ 5,{M,ff), (41) 

where we take ~ 24 to be the typical gauge coupling near unification. The 
unknown parameter 6^ should be of order one. Estimating Me// using a simple SU (5) 
model, we find a typical effective unification scale of 1 — 5xlO^''Ge^. This is only 
intended to a very rough approximation since a naively simple SU{b) model was used. 
Nevertheless, more complicated and realistic GUT models yield a unification scale in 
the same ballpark. It is generally true that our effective unification scale is about an 
order of magnitude or more greater than what is typically called the unification scale 
(~2xl0^6GeV). 

It may seem that our definition of an effective unification scale is rather artificial. 
However, it may be physically motivated by the following considerations. If indeed the 
standard model is embedded in some unified theory of gravity and gauge forces, then 
there may exist a phase at energies below the Planck scale which consists of a simple 
Lie group containing the supersymmetric standard model. In the absence of any 
gravitational corrections, the running gauge couplings certainly unify asymptotically, 
as this is the only case in which the higher group symmetry is fully realized up to 
arbitrarily high energies. Hence, the running couplings should only deviate from 
asymptotic unification by the gravitational corrections parameterized above. So, by 
the above reasoning, the effective unification scale should roughly correspond to the 
physical unification scale when the full (quantum gravitational) theory is considered. 

These results may have consequences for the paradigm of string unification. In 
particular, one problem of string unification |2Zj is that the couplings seem to unify 
at a scale [Mq^ ^ 2xlO^^GeV) about twenty times lower than the scale predicted 
by four dimensional heterotic string models (M^*^*""^ «i SxlO^^GeV). In the approach 
presented here, heavy threshold effects seem to push the effective asymptotic unifica- 
tion scale to roughly M^*""^. Despite the apparent success, this coincidence cannot 
be taken seriously until several questions are addressed in regards to this so-called 
string gauge coupling problem. First, the calculation of Mq""^^ j2H| was performed in 
the DR scheme, with the field theory step-function treatment generalized to strings. 
An analogous calculation for physical renormalization schemes is lacking, so it is 
difficult to compare our results with string predictions. Secondly, the asymptotic 
unification boundary conditions are probably not valid for many string models, and 
so the unification scale will be further changed by more general boundary conditions, 
as discussed in the next section. See for a recent string calculation of threshold 
corrections to grand unification. 
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5.4 More General Boundary Conditions 

The discussion of this section concerns the next-simplest boundary conditions after 
asymptotic unification. In particular, we will impose ai{Mu) = a2{Mu) = a^i^Mu) 
at scale Mu^Mh, for some h & H. As discussed in the previous section, one might 
expect Mjj to roughly correspond to the asymptotic unification scale, which we found 
to be roughly 5xlO^^GeV. However, we will consider Mu as an input and find the 
corrections for the strong coupling and gauge boson mass predictions. 

Before proceeding, there is a subtle point that should be addressed. Notice that in 
the previous section, we assumed that unification would have occurred asymptotically 
were it not for gravitational corrections. Hence, starting with a finite unification scale 
and then neglecting gravitational corrections, as we do in this section, does not seem 
logically consistent with what was done in the previous section. This observation 
is entirely correct, but the point is that indeed we are considering two orthogonal 
scenarios, one where a finite unification scale is obtained from gravity, and another 
where finite unification scale is obtained from non-trivial threshold corrections. The 
latter case may have its origin in stringy or gravitational physics, but nevertheless 
becomes manifest through purely field theoretic mechanisms. 

The corrections from imposing finite unification scale are straightforward to derive 
and can be stated in terms of the Aj = Af' + + 6i which we defined earlier. This 
gains an additional contribution and can now be written 

A, = Af + 5f + e,- Af , (42) 

where 




which is of exactly the same form expected from Eq. p5|) '^'^. Evidently, these are finite 
heavy threshold corrections in addition to the corrections from the heavy threshold 
splittings. Hence, the defined earlier will get an additional contribution from the 
Af^'s and is now 




Experimentally, — 1, as seen in Figs. |4I5] for typical gluino to wino mass ratios; 
this value can be easily adjusted for nonstandard sparticle spectra. This is our final 
formula which may be used to assess the experimental validity of gauge coupling 

^•^This is not obvious; one must work through the derivation to see that indeed the expected 
Lp — rjp correction terms do arise. 
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unification in any specific GUT model where the gauge group, superheavy mass ratios, 
and hght SUSY masses are given. 

Let us now consider the numerical size of these new threshold corrections. From 
Eq.(jni) one finds that 



which can be larger than heavy splitting corrections log jp- for values of Mu that are 
not too large. Hence, such corrections cannot be neglected. 

The value of M(/ is not fixed a priori, and corresponds to the physically meaningful 
energy where the couplings become equal due to the new nontrivial heavy threshold 
corrections. This complicates the analysis of unification by introducing another pa- 
rameter beyond those that are usually needed. However, this is to be expected, since 
a new physical phenomena (corrections arising from the virtuality of very massive 
particles) has been included. 

6 Conclusions 

We have developed a new way of looking at detailed predictions of gauge coupling 
unification which is more physically motivated than conventional approaches. In 
addition to a dramatic paradigmatic improvement, novel heavy and light threshold 
corrections are obtained, and the resulting corrections to unification predictions are 
presented for a general GUT model. A natural extension of this work is a thorough 
analysis and classification of various unified theories. By calculating the Bh constants 
and the heavy spectrum, one may exclude or verify the gauge unification of a given 
model. 



We wish to thank Ratin Akhoury and Eduardo de Rafael for reading a draft of 
this paper and providing useful suggestions. 



Here we will give explicit formulae for the pinch technique effective couplings reg- 
ularized using dimensional reduction (DRED), which will be denoted with a tilde. 
These effective charges will be similar to those constructed in for QCD, and 

in Jlljn] for the electroweak sector. However, we will extend these to the mini- 
mal supersymmetric case, which involves explicitly including another Higgs doublet. 




(45) 
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Appendix A : The Pinch Technique 
Self-Energy-Like Effective Charges 
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and regulating the loop integrals with dimensional reduction (DRED), as opposed 
to dimensional regularization (DREG), which is used in most non-supersymmetric 
settings. It is well known that DRED preserves both supersymmetry and gauge sym- 
metry. Also, the effective charges presented in were in the on-shell subtraction 
scheme (Qo = 0), whereas here we will need the result for arbitrary renormalization 
scale. In the appropriate limits our results reduce to those given in [T3j and [TEj . 

The charges are constructed using the pinch-technique (PT), which allows one to 
extract the universal self-energy function in non-abelian gauge theories, thus leading 
to gauge invariant effective couplings which 

• contain explicit and complete mass-threshold behavior and 

• reproduce the conventional massless beta function in the limit where masses 
can be neglected. 

At one-loop, the spin 1/2 and spin PT self-energies are trivially just the usual 
transverse vacuum polarization graphs. Only the graph with a gauge boson loop 
needs to have the self-energy-like part projected, as described briefly in section III, 
and in more detail in the references JHUHl- In calculating the following, we used both 
the direct diagrammatic pinch technique algorithm^H] and the dispersive derivation 
from physical cross sections [TB]. 

The PT effective charges naturally measure the self-energy-like propagation of a 
gauge boson and hence can be interpreted as measuring the real force between two 
fermions of arbitrary mass, analogous to the QED effective charge. The PT charge 
includes finite mass recoil effects that are missed in the heavy quark effective charge 
(the V-scheme). In fact, one may obtain the heavy quark potential in the appropriate 
kinematical limit of the pinch technique effective charge jl5j. The difference between 
the two are due to finite mass test-charge effects that are not present in the (V) charge 
but are in the (PT) charge. Different extensions of the PT effective charge beyond 
one-loop have been put forth [IZj [TH] [IHl , although it seems that the approach of [12] 
most closely matches the philosophy used here. A multi-loop generalization of this 
algorithm remains to be constructed. 

QCD Effective Charges 

The PT self-energy function for supersymmetric QCD, lis, can be used to define 
the effective coupling for supersymmetric QCD by 

53(g) = . (46) 

l + n3(Q,Qo) 

The function lis can be written down straightforwardly using Eqs. |5p9j . and the 
unsubtracted result is given by 



^N,{\og% - Cuv - t'ilii) 



3 ^ /i' 
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- -N,(L._{Q/m~) + log (m|/^2) _ _ 5/3 

- T.\{h (Q/m,) + log (mj//i2) - Cuv - 5/3 



E^(^o(Q/m~)+log(mf//.^ 
9 



~ C'l/v — 8/3 



(47) 



The four terms correspond respectively to the gluons, gluinos (^), Dirac quarks 
(g), and to complex squark doublets(g). For the scalars we will take the left and right 
components to be degenerate in mass since such complications do not change the 
unification predictions to any numerical significance. In any case, one may trivially 
treat the two separately. 

To relate the resulting effective charge to other schemes or observables one needs 
to use Eq.Q. 

Eq. ()47j) can be written in a more useful once subtracted form by relating the 
effective charges at different scales, leading to an expression governing the running of 
the charge given by 



n3(Q,Qo) = 



n3(Q) - n3(Qo) 

a3(<5o 



47r 



- Li{Qo/m~ 



o (^i (Q/ ^g) - L\ (Qo/ rrig] 

q -J 



1 



Lo(Q/m~)-Lo(go/ 



(48) 



Though the gluon contribution in Eq. (j47p looks simple, it is actually the most 
difficult piece to compute. As discussed in the pinch technique self-energy that 
contributes to the effective charge is gauge and scale independent, and indeed repro- 
duces the pure gauge term of the (3 function coefficient. This of course is not the case 
for the full pure gluon vacuum polarization, which is gauge and scale dependent and 
does not reproduce the correct (3 function. The non-trivial and important part of the 
result is the constant 64/33. This constant may be obtained using the pinch tech- 
nique and DRED. To translate to DREG one just subtracts 1/11 (from the so-called 
epsilon ghosts) to get the constant 67/33. For comparison, the heavy quark potential 
effective charge, ay, replaces this constant with 28/33 and 31/33 when using DRED 
and DREG, respectively. Consequently, the V-scheme doesn't satisfy the decoupling 
criterion of Eq.fllUj). This is just a reflection of the fact that infinitely heavy external 
quarks are used in the V-scheme calculation, thus rendering meaningless the limit 
where internal virtual gauge bosons acquire very large mass. 
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Notice that in the appropriate hmit the above reduces to the standard RG /3 
function coefficient, 

i^^on3 = 7;(9-^/)log^- (49) 

rrii^O 471" (^Q 

The Electroweak Sector 

The effective QED charge and the effective weak-mixing angle are obtained by 
diagonahzing the electroweak neutral currents and are given by 



and 



5(Q) = —59^ , (50) 

i + n^^(g,Qo) 



sUQ) = sim [i + ^i^i.ff'^;^"^, ) • (51) 

V s^(<5o) 1 + n^^(Q,Qo)/ 



For the matter sector, we will write only the subtracted PT self-energies, as it 
is now clear how to translate between schemes using the rip constants as described 
before. The quarks (g) and leptons (/), along with their scalar superpartners (g, /), 
yield 



Il-yy (matter) 



a{Qo) 



Y.\N,el(L.{Q/m,) - L.{Q,/m, 

q <J 



An 

I- H 

- E \N4{L,{Q/m~) - Lo(goK 

- Y.\(L.{Q/mi)-^{Q,/mi) 



2 



3 



(52) 



The electric charge of a particle p is denoted Cp. The analogous contribution of 
individual Dirac mass eigenstate matter fields to the 7Z self energy are given by the 
relation 

ng=(T^-^^)-VnS, (53) 



VAlp I ""JsC 

where p denotes any of the fermions or scalars above. 

The contribution of the charged vector bosons to the self-energies is more compli- 
cated than the matter multiplets. Similar to the QCD case, the non-abelian nature 
of the theory implies that W~^W~ loops (along with possible gauge dependent ghost 
and Goldstone boson loops), do not yield a gauge invariant result, and do not give the 
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appropriate contribution to the electroweak beta functions. The proper treatment in- 
volves calculating the self-energy like part of the one-loop e+e~^e+e" amplitude (or 
using any other fermions due to universality), including vertex and box corrections 
involving neutrinos. These contribute pinched parts which make the self-energy-like 
part gauge invariant and transverse. This calculation was first performed in ^3] , and 
then with dispersion relations in ^E], for uh = ^ Higgs doublets and renormalized in 
the on-shell scheme at Qo = 0. Here we need to extend these results to arbitrary uh 
and Qo, and would like to have the finite constants in the unrenormalized expression, 
including constant terms arising from using DRED instead of DREG. The most ef- 
ficient way to do this is to use the Feynman gauge ^ = 1, where bosons, — 
Goldstone bosons, and ghosts all propagate with —ig^u/ip"^ — ^w)- Hence, the 
only factors of transverse momenta arise from the three boson vertex, and so the box 
graph and several of the vertex graphs may be neglected, as they do not have pinched 
parts. Here, the dependence on the Higgs doublets comes only from the unphysical 
charged Goldstone scalars for the nu = 1 case, and also on charged Higgs for uh > 1. 
The result for the SU{2)l®U{1)y electroweak theory is 



l[iW,H) 



77 



q(Qo) 

47r 



2^(Li(Q/M^) + log {Ml/^l^) - Cuv - 64/33 



3 

+ E ( - o) (^o(Q/Mj + log (M,V/i') - Cuv - 8/3 

a=l 



(54) 



where the constant 64/33 is the same as appeared for the gluon self energy. The sum 
in the second line will be over mass eigenstate charged Higgs scalars; there will be 
one of these for each Higgs doublet in the theory. The first scalar (a = 1 in the sum) 
is an unphysical Goldstone boson that is eaten by the W^, and hence one identifies 
its mass to be Mi = Mw (in the Feynman gauge). The second charged scalar (a=2) 
is conventionally denoted by if^ in the MSSM, with mass M2 = Mh±. Additional 
Higgs doublets beyond the MSSM are not considered here so we can take tlh = 2. 
The function 

Ll{Q/m) = ^ (/5tanh- V"' (12 - /?') + - l) , (55) 



with (3 = ^1 + comes from the W+W' and ghost loops, the W+G' + W'G-^ 
loops, and the pinched self-energy-like part of the jWW vertex where the internal 
neutrino line is pinched. The Lq comes from the charged Goldstones and Higgs 
loops. As might be anticipated from the fermions and scalars, where, for example, 
lim^^oo -^i(Q/'"^) = 5/3 is the same constant as appears in the self-energy, we also 
have the nice property that 



M 



lim Li{Q/Mw) = 64/33. (56) 



Notice that in DREG this does not cancel the constant, which in that case is 67/33. 
With DRED regularization, all massive particles decouple, modulo divergent pieces. 
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from the unsubtracted self-energy-like expression. 

Letting the W-bosons eat the Goldstones by performing simple algebra in Eg. ()54|l . 
one finds the result written in terms of physical degrees of freedom, 



fliw,H) 



77 



Air 

+ E 

a=2 



7[Li{Q/Mw) + log (M^//x2) _ Cuv - 40/21) 



LoiQ/Ma) + log (M„7/i^) - Cuv - 8/3 



(57) 



The contribution of the physical massive gauge boson is characterized by 7Li 
(22/3)17+ (-l/3)Ls, explicitly given by 



L,{Q/m) = 2(3tanh-'l3-\l - (p' - l)/7) + (2/7)(/3' - !)• 



As expected, 



lim Lw{Q/m) = log 



lim Lw{Q/m) =40/21. 



(5J 



(59) 



Notice that Eq. (j3S|) precisely corresponds to Eq.(j2SI) for s = 1. 

The separation of pure gauge effects and those arising in the broken phase of the 
theory is useful, and allows us to immediately write down the analogous result for jZ 
without further calculation : 



{W,H) 



oj 



2^cI(l^{QIMw) + log (M^//i2) - Cuv - 64/33) 



1 1 

+ E ( - o (4 - ^)) (^o(Q/M,) + log (M.V/i^) - C^v - 8 



a=l 



(60) 



Finally, the wino and charged Higgsino, whose mixing is neglected, contribute 



n, 



H,W 
77 



and 



n 



H,W 



I \ 4 / 

(Q/m^) - Li_{QQ/mjj) j - -[Li{Q/m^) - LiiQ^/m^ 



4 / 

-i;cl[Li{Q/m^)-Li{QQ/m^\ 



(61) 



--(c^-l/2)(L.(g/m^)-L.(Qo/m^) 



(62) 



The SU{2)i®U{1)y effective couplings constructed from the above results are 

5^(Qo) 



l + n,(Q,Qo)^ 
33 



(63) 



for i = 1,2 and ai{Q) given in Eq.(j22I)- The PT self-energies are related by 



Hi — n^^ U^z 



C 



U2 = n^^ + — n^z. (64) 

Notice that fl^-^ = J2p n^^^ and Il^z = J2p n^z (as well as flj ) have the correct beta 
function coefficients, which are summarized below, and smoothly interpolate between 
all mass thresholds. The full mass-dependent beta functions may be obtained by 
differentiating the above expressions, but we will just give the massless limits, in 
order to make clear our conventions. 

The one-loop beta function coefficients are defined by the relations 

dsl dU^z a 

- ^-c^TfTTT^ = TZ(^^Z (66) 



c/logQ2 "'rflogg2 

?i = ■^(c^/377 - fi-^z) 

P2 = sl(3^^ + (3^z. (68) 



Pi = liclP,, - P,z) (67) 



One finds 



16 

/377 = ^^9 + Q-nH 



16 1 
P^z = -2Ng + -slNg + 6cl-nH{--sl), (69) 

thus leading to the MSSM beta function coefficients 

(pl) = (6]-nJ2] -nni 1/2 V (70) 




The two loop effects contribute to the running of the couplings through the terms 
9i{Q,Qo) in Eq.(jl7j). Since we do not have the full mass-dependent contributions, 
we will have to settle with the using the usual massless limits. These are explicitly 
determined by solving the two-loop renormalization group equations and are given 
by 

0^iQ, Qo) = E ^ log (1 + «,(go)§ log (QVQl)) , (71) 



where the beta matrix is 




^MSSM^_i 25 24 ). (72) 
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